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found to be significant and have to be taken into account for a proper prediction of the decay 
widths. 
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1 Introduction 


The announcement of the discovery of a new boson by the LHC experiments ATLAS [T] and 
CMS [23 has marked a milestone for particle physics. While the properties of this particle are 
consistent with the SM predictions, the uncertainties in the experimental data still leave enough 
room for interpretations in extensions beyond the SM. Among these, supersymmetric (SUSY) 
models [3] certainly belong to the best motivated and most intensely studied ones. In particular 
the Next-to-Minimal Supersymmetric Extension (NMSSM) [4] provides with the introduction of 
an additional complex superfield S a dynamical solution to the n problem [5] when the singlet 
field acquires a non-vanishing vacuum expectation value. Because of new contributions to the 
quartic coupling A, with which S couples to the Higgs doublet superfields H u and H c i-> the 
tree-level mass value of the lighter MSSM-like Higgs boson is enhanced. In consequence less 
important radiative corrections are required to shift the mass value to the measured value of 
125 GeV and therefore smaller stop masses and/or mixing are necessary, so that the fine-tuning 
is reduced m- 

After electroweak symmetry breaking (EWSB) the Higgs sector of the NMSSM consists of 
seven physical Higgs bosons. In the CP-conserving case, which we assume to be valid here, these 
are three neutral CP-even, two neutral CP-odd and two charged Higgs bosons. The discovery 
of all Higgs particles is challenging though not impossible at the high-energy option of the 
LHC. In [8] we investigated the discovery prospects for the NMSSM Higgs bosons during the 
13 TeV run of the LHC and gave benchmark scenarios that feature Higgs-to-Higgs decaysQ If 
kinematically allowed also decays into supersymmetric particles can become important, as is well 
known for the Minimal Supersymmetric Extension of the SM (MSSM) [37/38;]. The one-loop 
SUSY-QCD corrections to the decays into stops and sbottoms of the MSSM Higgs bosons have 
been calculated in [39H4T] and can change the decay widths by more than 50%, especially near 
the threshold. The SUSY-QCD corrections have been reanalyzed in [32] • The full electroweak 
(EW) one-loop corrections to the pseudoscalar decays into squarks have been provided in [33] 
and have turned out to be significant. Equally, the decays of heavy squarks into lighter ones 
and a Higgs boson can dominate in a wide range of the MSSM parameter space due to the large 
Yukawa couplings and stop and sbottom mixing [33]. The SUSY-QCD corrections at next-to- 
leading order (NLO) are of the order of a few ten percent [35] and mostly negative. The one-loop 
EW corrections to the decays with a pseudoscalar in the final state are significant [43] . The full 
one-loop corrections for the complex MSSM have been discussed in [46] . 

The proper interpretation of the experimental data and, once SUSY has been discovered, the 
aim to pin down the underlying model and distinguish e.g. the NMSSM from the MSSM, require 
precise predictions both for the parameters of the model and for the observables like e.g. NMSSM 
Higgs boson production and decay rates m- The higher order corrections to the CP-conserving 
NMSSM Higgs boson masses and self-couplings have been given in [48H58] and m, respectively. 
Additionally, there are several codes available for the evaluation of the NMSSM mass spectrum 
from a user-defined input at a user-defined scale, like NMSSMTools [59H61] which can be interfaced 
with SOFTSUSY [621163] . the interface of SARAH [571I64W6T] with SPheno [681169] . and finally SARAH 
which has been interfaced with the recently published package FlexibleSUSY [701171] . Recently, 
NMSSMTools has been extended to include also the CP-violating NMSSM [72] . In our Fortran 
package NMSSMCALC m we have included in the CP-conserving and CP-violating NMSSM the 
full one-loop and the order 0(a.ta s ) corrections to the NMSSM Higgs boson masses [561741175] and 

x For other recent studies on NMSSM Higgs boson phenomenology, see Refs. miTTl h>(il . 
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the state-of-the-art higher order corrections to the decays. These include in the CP-conserving 
case Higgs decays into stops and sbottoms the SUSY-QCD corrections of [52] which have been 
adapted from the MSSM to the NMSSM case. Very recently, neutral Higgs production through 
gluon fusion and bottom-quark annihilation including higher order corrections has been discussed 

hr [76] . 

With this work we take another step in improving the predictions for the NMSSM Higgs 
sector. We provide both the NLO SUSY-QCD and the full one-loop EW corrections to the 
decays of a pseudoscalar NMSSM Higgs boson into stops as well as to the decays of the heavier 
stop into the lighter one and a pseudoscalar. 

The paper is organized as follows. In section [2] we introduce the NMSSM Higgs sector and 
set our notation. The tree-level decay of a pseudoscalar Higgs into stops is discussed in section 
El before we present in section El the order 0(a s ) SUSY-QCD corrections and in section El the 
one-loop EW corrections. The numerical analysis, including the discussion of heavier stop decays 
into a pseudoscalar and lighter stop final state, is performed in section [G] We conclude in section 

m 

2 The NMSSM Higgs Sector 

The NMSSM Higgs potential is obtained from the NMSSM superpotential, which we assume to 
be scale invariant, the soft SUSY breaking terms and the D-term contributions, which are the 
same as in the MSSM. In terms of the superfields H u and H d , coupling to the up- and down-type 
quarks, respectively, and the singlet superfield S the superpotential reads 

Wnmssm = Wmssm — + -K.S 3 , (2.1) 

where i. j = 1,2 are the SU(2)l indices and we have introduced the totally antisymmetric tensor 
eij with ei 2 = 1. Working in the CP-invariant NMSSM, the dimensionless parameters A and k 
are chosen to be real. The MSSM superpotential Wmssm is given by 

Wmssm = e# [y e H l d VE c + y d H l d QW c - y u H i u Q j U c } , (2.2) 

with the quark and lepton superfields and their charge conjugates, indicated by the superscript 
c, denoted by Q,U C , D c , L an E c . The color and generation indices have been suppressed in 
Eq. (12.21) . We neglect generation mixing of the quarks, so that the phases of the Yukawa 
couplings y d , y u and y e , which in general are complex, can be reabsorbed by a redefinition of 
the quark fields. The soft SUSY breaking NMSSM Lagrangian involving the Higgs doublet and 
singlet component fields H u , H d and S reads 

£-soft = £-soft, mssm ~ "4 , |5'| 2 + ( e ij^A\SH d H-[ — — k,A k S 3 + h.c .) , (2-3) 

with the soft SUSY breaking MSSM Lagrangian 

C-soft, mssm = -m 2 Hd H\H d - m 2 Hu HlH u - rri^Q^Q - m\L ] L - m 2 ^ R U* R u R - rn 2 fji d* R d R 
-m% R e* R e R - (e i: [y e A E H^d* R + y d A D H^d* n - y u A v u* R ] + h.c.) 

— + M 2 W k W k + M 3 GG + h.c.) . (2.4) 
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The gaugino fields are denoted by B , Wk (.k = 1,2,3) and G, and the left-handed squarks and 
sleptons are arranged in doublets denoted by Q = ( UL,di ) T , L = while the right- 

handed fields are denoted by ur, d r and e.R. The soft SUSY breaking mass parameters 

of the scalar fields X = S, H d , H U ,Q,U , D, L, E are real, while the gaugino mass parameters 
Mi, M 2 and M 3 and the soft SUSY breaking trilinear couplings Ay (Y = A ,k,U,D,E) are in 
general complex. In the CP-conserving case assumed here, they are, however, real. Again, the 
respective quark and lepton superfields are understood to refer to all three fermion generations. 
Note that we have set soft SUSY breaking terms linear and quadratic in the singlet field S to 
zero. 

After expanding the Higgs fields about their vacuum expectation values (VEVs) v u , v d and 
v s , chosen to be real and positive, 


H d = 


(v d + h d + ia d )/V 2 \ 

hi ) 


H u = 


K 

u Y h u Y i&u)/V% 


S = 


Y h s Y id s 

~V2 


(2.5) 


the Higgs mass matrices for the three scalar, two pseudoscalar and the charged Higgs bosons 
can be derived from the tree-level scalar potential. The mass matrix decomposes into two mass 
matrices for the CP-even and the CP-odd Higgs fields. The squared 3x3 mass matrix Mg for 
the CP-even Higgs fields can be diagonalized through a rotation matrix 1Z S which yields the 
CP-even mass eigenstates Hj (i = 1,2,3) as 

(H U H 2 ,H 3 ) T = 'Jl s (h d ,h u ,h s ) T . (2.6) 


The Hi are ordered by ascending mass, Mh 1 < M # 2 < M# 3 . In order to obtain the CP-odd 
mass eigenstates A\, A 2 and the massless Goldstone boson G first a rotation VP to separate G 
is applied, and then a rotation 1Z P to obtain the mass eigenstates 

(A 1; A 2 , G) t = TZ p (a , a Sl G) T = TZ p U G (a d , a u , a s ) T , (2.7) 


which are ordered such that Ma 1 < M 4 2 . 

The minimisation of the Higgs potential V requires the terms linear in the Higgs fields to 
vanish in the vacuum. The corresponding coefficients, which are called tadpoles, therefore have 
to be zero. The tadpole conditions for the CP-even fields can be exploited to replace m 2 Hu , rn 2 Hd 
and m 2 s by the tadpole parameters th d ,th u and th s . Replacing the SU(2)l and U(l)y gauge 
couplings g and g' and the VEVs v u and v d by the electric charge e, the gauge boson masses 
M\y , Mz and by tan/3, the tree-level NMSSM Higgs sector can then be parameterized by the 
twelve parameters 

thu, th d , t hs , e, Myy, Mf, A, At, Ax, A K , tan/3 = (H°)/(H%) and g eS = X(S) . (2.8) 

The VEVs of the neutral components of the Higgs fields are denoted by the brackets around 
the corresponding fields. The sign conventions for A and tan /3 are chosen such that they are 
positive. The At, A \, A K and // e ff on the other hand can have both signs. Note also, that the 
parameter A\ can be traded for the charged Higgs boson mass M H ±, which we will do in the 
following. From now on we will drop the subscript ’eff’. Note that the inclusion of higher order 
corrections requires also the soft SUSY breaking mass terms for the scalars and the gauginos as 
well as the trilinear soft SUSY breaking couplings. 
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3 The Tree-level Decay Width 


We start by discussing the tree-level decay width of a pseudoscalar Higgs boson At (i = 1,2) 
into a pair of stops. The stop mass matrix in the interaction basis reads 


M} = 


m 

m 


LL 

2 

RL 


m 

m 


LR 

2 

RR 


(3.9) 


with 


12 

m LL = m 1j + m t + M | cos 2/31- — - sin 2 9\y 


2 3 

2 , 


m RR = + mf + cos 2/3 sin 0m 

t-R 3 


- n cot /3) , 


(3.10) 

(3.11) 

(3.12) 


in terms of the soft SUSY breaking mass parameters tUq and the soft SUSY breaking 

trilinear coupling At, the higgsino mixing parameter /r, the top and the Z boson masses mt and 
Mz, the mixing angle f3 and the Weinberg angle 9\y. The /i parameter is generated dynamically 
in the NMSSM and given by 


V = 


A v s 
sft 


(3.13) 


The stop mass matrix is diagonalized by 


TV = 


cos 6i sin 9i 
— sin 9i cos 9 1 


(3.14) 


yielding the stop mass eigenstates U (i = 1,2) as 


ti = tVj, 




(3.15) 


where s = L,R and for the squark masses we have < m^ 2 by convention. The mixing angle 
9i and the squark masses are given by 


tan 9i = 


^ m LR 


m LL - m RR - \/( m LL - m RRY + Am LR 


(3.16) 


and 


m~ = - 

h ,2 2 


+ min -F \/ ( m l L ~ m RR ) 2 + Am LR 


(3.17) 


In case of the pseudoscalar only the coupling to two different stop mass eigenstates is non¬ 
vanishing. For the tree-level decay width r LO we have 


T lo (T ->■ Ut 2 ) = 


3A 1 / 2 (M\., m~ , m~ ) 


87 rM], 


£ z -i° 

3 = 1 


12 


(3.18) 
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where X(x, y,z ) = (x — y — z) 2 — 4 yz is the two-body phase space function and the coupling G^ 2 
(j = 1. 2) of the pseudoscalar Aj to the stops reads 


r u _ am 

Aj 2 M w 


A t 

tan /3 



Kfi + 


Xv 1 
y/2 tan /3 



(3.19) 


where v is the VEV given by v 2 = v 2 + v ^ and 1Z p k are the elements of the rotation matrix 
defined in Eq. (12.711 . In particular, they are the tree-level mixing matrix elements. In the 
kinematics of the decay, however, i.e. for the external Higgs field, we use the two-loop corrected 
Higgs boson masses at order O(cttot s ), which include the full EW corrections at one-loop order. 
The renormalization of the Higgs fields and the computation of the mass corrections have been 
described in Refs. usmsizs]. We follow the conventions of these papers, to which we refer 
the reader for more details. In order to ensure the on-shell properties of the external Higgs 
field, which in the calculation of the Higgs mass corrections has been renormalized in the mixed 
on-shell DR scheme, the finite wave function renormalization factors Z jTT] have to be taken 
into account. The application of the factor Z to the tree-level matrix 1Z P (in G 1 ^.) leads to the 
rotation matrix 1Z P ' 1 which, modulo the Goldstone boson G , rotates the interaction eigenstates 
a and a s to the loop corrected mass eigenstates A\ and A^. cf. [56], 


— ZyK?i , i,j — Ai,A%,G , l — a,a s ,G . (3.20) 

These and the loop-corrected masses are taken from the Fortran code NMSSMCALC 1731 . in which 
we choose the on-shell (OS) renormalization for the top/stop sector in order to be in accordance 
with the renormalization scheme chosen later on both in the SUSY-QCD corrections and in the 
electroweak corrections. 

In Eq. (13.18]) we have summed over both possible final states tit^ and £*0- In the MSSM, 
i.e. leaving out the singlet contribution oc A in the coupling G^ 2 , the decay width is proportional 
to m 2 (y + At/ tan /3) 2 /M^. For small values of tan/3 and not too heavy pseudoscalars, the 
decay into stops can compete with and even dominate over the decays into top quarks and 
into charginos and neutralinos. In the NMSSM, this statement has to be taken with caution, 
however, as the singlet component in the coupling G^ 2 , depending on the scenario, can come 
with both signs and hence increase or decrease the decay width. 


4 SUSY-QCD Corrections 

The SUSY-QCD corrections for the NMSSM pseudoscalar A, decay width differ from the ones 
of the MSSM (39l |41~1 solely in the tree-level coupling to the stops G^ 2 . We shortly repeat them 
here for completeness and in order to introduce our renormalization scheme. 

The virtual corrections at order 0(a s ) to the pseudoscalar Higgs decays into stops consist of 
loop diagrams with a gluon, respectively, gluino exchanged in the A{t\t 2 vertex and of a contri¬ 
bution involving the four-squark vertex, cf. Fig. [T] (upper). Note that the mixing contributions 
due to off-diagonal self-energies are absent in the case of pseudoscalar Higgs decays, as At only 
couples to different stop mass eigenstates. The computation of the virtual diagrams leads to 
ultraviolet (UV) and infrared (IR) divergences. We work in dimensional reduction [78] , which 
preserves SUSY at the one-loop level. The fields and couplings are then treated in 4 dimensions 
while the loop integral is performed in D = 4 — 2e dimensions. The UV divergences appear as 
poles in e and are canceled by the wave function counterterms and the counterterm renormal¬ 
izing the Ait\t 2 interaction. The infrared divergences are regularized by the introduction of a 
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Figure 1: Diagrams contributing at NLO SUSY-QCD to the virtual (upper row) and real corrections (lower row) 
of the decay F(A —► txtz)- 


fictitious gluon mass The IR divergences left over after renormalization are canceled after 
adding the real corrections. These consist of the radiation of an additional gluon off the final 
state stops and are shown in Fig. Q] (lower). 

The one-loop corrected decay amplitude can be written as 


with 


-pNLO _ 'pLO I -p(l) 

1 QCD - 1 + 1 QCD > 


r 


(i) 

QCD 


= Re 




47rM^. 



q 12 * 


v x> 


A 


QCD 


\k = l 


(4.21) 


(4.22) 


where again the Z factors appear to ensure the on-shell properties of the external loop-corrected 
Higgs field. The A^ D are given by the sum of the virtual, real and counterterm contributions 
A^ fc and Arespectively, 

A^° = A^ + A^ + A^ . (4.23) 

From now on a factor Cfoi s /{A it) with Cf = 4/3 is factorized out and already included in 
Eq. (14.221) . The virtual corrections receive contributions A^ from the gluon exchange diagram, 

A^ from the gluino exchange diagram and A 44 from the diagram involving the 4-squark vertex. 
They are given by 


A 




+2(m| + m| - M\ k )C 0 {ml , M \, m|; C, m h , m h ) 


(4.24) 


for the gluon exchange, 


= gmtKki 
Ak Mw tan P 


(m-g - m t sin(2%)) -B 0 (m~ ; m g , m t ) + ( m- g + m t sin(2%)) R 0 (m|; m g ,m t ) 


- (sin(2%)m t (m t - i - ) + m g M A \ C Q {m i ,M Ak ,m i ]mg,m t ,m t ) 


(4.25) 
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for the gluino exchange, and 



G A k 


(4.26) 


for the 4-squark vertex diagram, where Bq and Cq are the Passarino-Veltman scalar two- and 
three-point functions m, cf. App.[A]for their definitions. 

The counterterm corrections consist of the renormalization of the external squark wave func¬ 
tions Z^ i (j = 1,2) and the renormalization of the Akt\t 2 interaction vertex. Note that the 
wave function renormalization of the pseudoscalar does not contribute at order 0(a s ) and hence 
to the SUSY QCD corrections. The parameters A, v, v s , tan/3, M\y, Mz and e are not renor¬ 
malized by the strong interaction, so that the counterterm reads 


A ct 
A A k 


12 


= ~Ga 
2 A 


m,i,+Mi,h 


) + 


dmt 


Suit + 


d ^L 

dA t 


8 At 


(4.27) 


leading to 


= ^l(5Z ilh + SZ, 


5m t 


fn) + G 12 U " H 


t 2 t 2 J 


T A k 


+ 


gm t 


v p 

'Del 


mt 2 M\v tan f3 


8A t 


(4.28) 


The counterterm 5At of the trilinear coupling is given by the quark and squark mass counter¬ 
terms, dmt and 8m^ , and the mixing angle counterterm 56 £, 


5 At = 


2 mt 


(ml-ml) 


2 cos(2 6i)56i — sin(2@j 


dmt 

) -- 

m t 


+ 2sm(20 t -)(m fl <ym fl - m^Sm^) 


.(4.29) 


In the (s)quark sector we adopt OS renormalization with the quark and squark masses defined as 
the poles of their respective propagators and the squark wave function renormalization constants 
defined such that the residues of the poles are equal to one. Defining the following structure for 
the quark self-energy, where Vl,r denote, respectively, the left- and right-chiral projector, 

Z t (p 2 ) = ^(p 2 )V L + jffif ( p 2 )Vr + m t Z Ls (p 2 )V L + m t Z Rs (p 2 )V R , (4.30) 


we have for the top quark mass counterterm 


1 , 


5m t = -Re(mtSf(m^) + miEf(m?) + Sf s (m?) + Sf s (m^)) 
The squark mass and wave function counterterms are given by (j = 1, 2) 


5m\ = ReSr,r,(m?) 


dT, i . i .(p 2 ) 
= -Re 33 


5Zr 


dp 2 


p 2 =m~ 


(4.31) 


(4.32) 

(4.33) 


In Eq. (14.3211 the denote the diagonal parts of the squark self-energies. The diagrams 

contributing at order 0(a s ) to the squark and quark self-energies are depicted in Fig. [2] (upper) 
and (lower), respectively. 
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Figure 2: Diagrams contributing at NLO SUSY-QCD to the squark self-energies (upper row) and to the quark 
self-energies (lower row). 


The mass counterterms read 
5m t 


= —Re 


m t 


m t 1 


mt, C) + Bi(mf ; m,g, m t -J + -Bi(m 2 ; mg, m t - 2 ) + 4:B 0 {mf;m t ,() 
m, 


+ sin(2%)— - £ 0 (mf; rrig, m t -J) 

xn t 


(4.34) 


= Re 


1 


2mtmj sin(20 ? )R o (m|;m t , m^) + - [(l + cos 2 (2%)) 

+ sin 2 (2%) Ao ( m t 2 )] “ 2m |-^o(m|jm^C) - A 0 (mg) - A 0 {m t ) 
+(m 2 - m 2 - ; m t ,nig) 


(4.35) 


m- t2 5m i2 


= Re 


- 2m t mgsin(2%)R 0 (m|;m t , nig) + -[(! + cos 2 (2%)) A 0 (m£ 2 ) 


+ sin 2 (2%)^ 0 (m|- 1 )] - 2m|R 0 (m|; m ?2 , C) - A 0 (mg) - A 0 (m t ) 

+(m? - m? - ml)B 0 (m~ ;mt,nig) 


(4.36) 


where B\ is the coefficient of the two-point tensor integral of rank one and Ao denotes the scalar 
one-point function, cf. App. [Aj The wave function corrections can be cast into the form 


sz t!t i = Re 


- 4 m t m~g sin(26* t -)Ro(m 2 ; m t , nig ) + 2(m~ + m 2 - m~ )B' 0 (m 2 r ; m t , nig) 


+2R 0 (m~ ; m ?1 , C) - 2 B 0 (m? ; m t , m- g ) + 4m 2 ^(m 2 ; m h , C) 


2 r' ., 


^- 2 t 2 = Re 


4m t mj sin(2%)Bo(m 2 ;m t , mg) + 2(m| + m 2 - m\ )B' Q (m 2 r -m t ,m-g) 


+2 B 0 (mr ; m fa , C) - 2B 0 (m t - ; m t , m § ) + B 0 (m f ; m t ~ 2 , C) 


(4.37) 


(4.38) 


Here B' 0 (k 2 ;m 2 , m 2 ) denotes the derivative with respect to /t 2 . The mixing angle counterterm 
is renormalized as 


1 


1 


<5% = - (SZ^) = —^^ 

2 1 ^ 2 (mf — m? 

v Cl 


Re i^ht2( m V + ^t 2 i 1 ( m V\ > 


(4.39) 


* 2 ' 

























where denotes the respective squark self-energies, so that 

1 


Mi = 


2 2 
m~ — m~ 

1 1 C2 


Re 


sin( 2 %) cos( 2 %) ^ 0 (m|) - A(m?)) 


+2m t mgCos(26» t -) f R 0 (m|; m t , m^) + 


(4.40) 


The real corrections finally in terms of dilogarithms read 


= 


2g A fc r 
avs L 


{M\ k - m? - m|) ( - 2 log A) log 


C M A k m h m ~ t2 


+ 2 log 2 /3 0 


- log 2 A - log 2 /3 2 + 2Li 2 (l - /3 q) - Li 2 (l - ^i) - Li 2 (l - /3|)) 
+ 2 A 1 / 2 log + 4A 1 / 2 + ( 2 M| fc + m? + m|) log A) 

+( M A t + 2m? ) log /3 2 + (Mi + 2m? ) log A 


(4.41) 


where in the two-body phase space function A (M^ fc , m 2 ^, m 2 ^) we have neglected the arguments 
for better readability. We have furthermore introduced 


A) = 


Pi = 


P 2 = 


Mi — m? — m? + A 1 / 2 


ti 


ti 


2m il m i2 


Mi — m? + m? — A 1 / 2 


ti 


^2 


2M Afc m t - 2 

M? + m? — m? — A 1 / 2 


(4.42) 
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5 The One-Loop Electroweak Corrections 

The NLO electroweak corrections consist of the virtual corrections to the vertex and the coun¬ 
terterm contributions to cancel the UV divergences. Besides the top and stop fields and the 
parameters specified below, again the Higgs field needs to be renormalized, as for the external 
Higgs field we use the two-loop corrected Higgs boson mass at order 0(ata s ), including the full 
EW corrections at one-loop order. In the loops, however, the tree-level masses for the Higgs 
bosons have to be used so that the UV divergences are canceled properly. 

For the NLO EW corrections, in addition the mixings 

-> kh) (5.43) 

of the decaying CP-odd Higgs boson A with the Z boson and the Goldstone boson G have to be 
included. The matrix element for the EW corrected pseudoscalar Higgs decay into stops hence 
reads 


M(Ai 


tit 2) — ^ Z 
3 = 1 


v 


-G%+A™) 


+ 5M ; 


G,Z 


mix ,i ’ 


(5.44) 
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Figure 3: Generic diagrams contributing to the electroweak corrections of the decay r(A; —> t\t 2 ) with i,j, k,p = 
1,2,1 = 1, 2, 3 and m, n = 1,5. 


where A EW represents the sum of the 1-particle irreducible (1PI) diagrams A^’ EW 
to the EW virtual corrections of the vertex and of the counterterms A^ T ’ EW , 

A EW = a v,ew + a ct,ew 


contributing 


(5.45) 


Due to massless photons in the loops we also encounter IR divergences. These are canceled by 
adding the real corrections A^ ,EW , where a photon is radiated off the final state stop lines. We 
hence have for the EW corrected decay amplitude 


~^ ^1^2) — 


3A 1 / 2 

87 


E Z « G 


3=1 


+ E Z « G $ - E z -* a 


R, EW 
A k 


< 3 =1 


\fc=l 


-2Re [ ■£ Z‘fi'1' (E 

<3=1 \k=l J 


(5.46) 


Again we have dropped the arguments in the two-body phase space function A. In the following 
we will discuss the individual contributions. The virtual corrections consist of the 1PI diagrams 
given by the triangle diagrams with scalars, fermions and gauge bosons in the loops, as shown in 
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*), {H ± /G ± , b ), {tj, x°J, (h, Xp) 


Figure 4: Diagrams contributing to the electroweak squark self-energies (first two rows) and quark self-energies 
(last row) with i,j,k,p= 1,2, l = 1, 2, 3 and m = 1,5. 


the first two rows of Fig. [3l and of the diagrams involving four-particle vertices, cf. Fig. [3] (last 
row). For better readability, for the scalars S appearing in the loops we only listed the particle 
types but not the combination of scalars that are allowed by the theory for the various vertices. 
Let us remark, however, that in the four-particle vertices the scalar-Higgs—Goldstone-boson—2- 
stops coupling Hi — G — i\ — and the scalar—pseudoscalar—2-stops coupling Hi — A p — i\ — 

{l = 1,2,3 , p = 1, 2) are new compared to the MSSM. The former is due to the singlet admixture 
in Hi, the latter is proportional to the NMSSM specific coupling A. The diagrams have been 
generated with the Mathematica package FeynArts 3.6 |8U] and evaluated with FormCalc 7.3 [ST] 
in two independent calculations. The integrals have been computed with LoopTools 2.7 [8 Tj . 
The results of both calculations agree and have been cross-checked against a third calculation, 
that did not use any of the tools to evaluate and simplify the amplitudes, and which takes the 
loop functions from HDECAY [821183] and SDECAY [83]. The UV divergences encountered in the 
computation of the virtual corrections are canceled by the counterterms that are the sum of the 
stop wave function corrections and of the counterterm renormalizing the Ajtit? interaction, 


a c T , ew = a ct,„ a ct> 

Si j Si j 


(5.47) 


Because of the antisymmetric coupling the stop wave function corrections are given by 


G 12 

A%’ W = —£ L (8Z flh +6Z hh ). (5.48) 

The stops are renormalized on-shell, with the renormalization conditions given in Eqs. (14.321) 
and (14.331) . The diagrams, that contribute to the here required electroweak self-energies are 
displayed in Fig. [3] (upper two rows). 
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For the one-loop EW corrections the vertex counterterm reads 

*ct,» = _£-d2 (^L + _ 5M w \ _ gm t 

Aj Aj V 9 m t M\y ) 2M\y 

5X 5v <5tan/3\ Xv P 
X + v tan/3 / y/2tanj3 j2 _ 




<5 tan/3 \ A t \ 
tan /3 J tan j3 J 

(5.49) 


The individual counterterms are derived from the renormalization of the input parameters. We 
follow Ref. [56] and apply the same renormalization scheme which mixes OS and DR conditions 
as defined there. For the vertex counterterms the relevant input parameters are the W and Z 
boson masses M\y and M z , the electric charge e, tan/3, A and The parameters that can be 
related to physical quantities are renormalized OS, the remaining ones DR. Together with the 
OS-renormalized top/stop sector, we have the following set of parameters to be renormalized, 


tan/3, A, v s ,M z ,M w ,e, ra t , , m f - 2 , % . (5.50) 

DR scheme OS scheme 

The coupling g and the VEV v appearing in Eq. (15.491) are given in terms of these parameters 
by 


eM z 

yJ M z ~ 


and 


_ 2 M\y j~ M\y 


(5.51) 


from which their counterterms can be derived. The details of the renormalization of the counter¬ 
terms for the first six input parameters can be found in [[5B], so that they are not repeated here. 
The formulae for the OS renormalization of the top and stop masses are given in Eqs. (14.3111 - 
(14.33)1 . The squark and quark self-energies for the EW one-loop corrections are depicted in 
Fig. [4[ At EW one-loop order the counterterm for At now reads 


SA t = 


9 


9 


6 tan /3 


+ 


tan/3 \ /r tan/3 J ' 2 mt - 
+2 sin(2(%) {m- tl dm ^ - m~ t2 5m~ t2 ) 


(m| - ml) ( 25%cos(26> i ~) - sin(2%) 




Smt 

m t 


(5.52) 


The counterterm for the mixing angle 6i is renormalized as in Eq. (I4.39p . however with the 
self-energies given by the diagrams shown in Fig. [I] (first two rows). 

The diagrams for the contributions to the electroweak corrections stemming from the mixings 
of the pseudoscalar At with the Z boson and the Goldstone boson, 5 M v , are shown in Fig. [5[ 
Using the Slavnov-Taylor identity [85], on can obtain the mixing contributions through 


= 


12 

M z 


S AiZ 




(5.53) 


where Z^AiZ denotes the renormalized At-Z mixing self-energy and Gj? is the Goldstone boson 
coupling to the stops, 


r ,i2 _ 9mt f _ 

G 2 M\y \ 1 tan /3 J 


(5.54) 


2 Additionally, for the renormalization of the Higgs sector in the computation of the higher order corrections to 
the Higgs boson masses, we have the tadpole parameters, the mass of the charged Higgs boson M H ±, the NMSSM 
parameter k and the trilinear coupling A K , which need to be renormalized at loop level, cf. [56 . 
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Figure 5: Generic one-loop diagrams contributing to the mixing SM^r A of a pseudoscalar Ai with the Z and 
Goldstone boson. 


Note, in particular that the external momenta have to be set to the tree-level mass in order 
to ensure to maintain gauge invariance. 

The last piece which is missing in the decay, in order to also get an IR finite result, is the 
real corrections term A^_ EW . This is the same as for the QCD corrections, but with the gluon 
replaced by the photon. In the formula for the QCD corrections, Eq. (14.410 . this means that 
the coupling and color factors have to be replaced accordingly, i.e. 


a r,ew 

^A k 


167T 2 


Ai 


(5.55) 


with A^ fe given in Eq. (14.411) . 

The full NLO decay width including the SUSY-QCD and -EW corrections is then given by 


pNLO, full 


-pNLO I -pfl) 

1 EW T 1 q CD , 


(5.56) 


where Tg^P, defined in Eq. (15.461) . includes the leading order decay width and Tqq D has been 
defined in Eq. (14.211) . 


6 Numerical Analysis 

For our numerical analysis we first perform a scan in the NMSSM parameter space in order to find 
scenarios that are in accordance with the LHC Higgs and SUSY data. The compatibility with the 
LHC Higgs data has been checked by using the programs HiggsBounds [86j and HiggsSignals 
[57]. The effective couplings of the NMSSM Higgs bosons, normalized to the corresponding SM 
values, as well as the masses, the widths and the branching ratios of the NMSSM Higgs bosons, 
which are required as inputs for these programs, have been obtained from the Fortran code 
NMSSMCALC m ■ The loop induced Higgs coupling to gluons normalized to the corresponding 
coupling of a SM Higgs boson with same mass is obtained by taking the ratio of the partial widths 
for the Higgs decays into gluons in the NMSSM and the SM, respectively. These include the QCD 
corrections up to next-to-next-to-next-to leading order in the limit of heavy quarks [8814971 and 
squarks [981(99]. taken over from the SM, respectively, MSSM case, while the EW are unknown 
for the SUSY case and hence consistently neglected also in the SM decay width. The stop 
mass values have been chosen such that they are not excluded by present ATLAS [10011102] and 
CMS [10311104] searches. The squark masses of the first two generation are heavy enough not to 
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be in conflict with LHC data. The SM input parameters that we use are [105 . 106 ' 


a(M z ) = 1/128.962, a^ 3 4 S (M z ) = 0.1184, M z = 91.1876 GeV, (6.57) 

M w = 80.385 GeV , rn t = 173.5 GeV , mp(mp) = 4.19 GeV . 

The running ct!? R used in NMSSMCALC is obtained by converting the that is evaluated with 

the SM renormalization group equations at two-loop order, to the DR scheme. The light quark 
masses, which have only a small influence on the loop results, have been set to 

m u = 2.5 MeV , = 4.95 MeV , m s = 100 MeV and m c = 1.42 GeV . (6.58) 

As renormalization scale we choose the SUSY scale M s , which we set 

M s = s/ m Q 3 m t R = VR ■ (6-59) 

6.1 Pseudoscalar Higgs Boson Decays into Stop Pairs 

In this subsection we present the impact of the SUSY-QCD and -EW corrections on the decay 
of a heavy pseudoscalar Higgs boson into a pair of stop quarks. The parameter point, which we 
have chosen from the set of parameter points that survive the LHC constraints, is given by the 
soft SUSY breaking masses and trilinear couplings 

m un,c R = m d R ,s R = m Qi, 2 = m Zi, 2 = = 3 TeV, m~ tR = 536.43 GeV, 

ttIq 3 = 594.61 GeV, mj )R = 1285 GeV, m^ 3 = 255.53 GeV, nrif R = 1499 GeV, 

A t = 1418 GeV, A U)C = 1435 GeV, A d>Smb = -66.68 GeV, A e ^ r = -91.76 GeV, 

Mi = 111.73 GeV, M 2 = 395.86 GeV, M 3 = 1370 GeV (6.60) 

and NMSSM specific input parameters 

A = 0.629 , k = 0.223 , A K = -543.53 GeV , /r eff = 452.61 GeV , 

tan/3 = 1.969, M H ± = 1024 GeV . (6.61) 

This results in the two-loop corrected mass Ma 2 of the heavy pseudoscalar A 2 and the stop 
masses rn^ 2 , 

Ma 2 = 1012 GeV , m h = 280.78 GeV and m t ~ 2 = 709.07 GeV . (6.62) 

We follow the SLHA format [107] , in which the parameters A, k. A k , // e ff, tan f3 as well as the soft 
SUSY breaking masses and trilinear couplings are understood as DR parameters at the scale 
fin = AfJl, whereas the charged Higgs mass is an OS parameter. As input for our computation we 
take the soft SUSY breaking trilinear stop coupling At, however, consistently as OS parameter. 
The conversion from the DR to the OS scheme is done within NMSSMCALC and yields 

A° s = 1435 GeV . (6.63) 

3 For tan/3 this is the case only, if it is read in from the block EXTPAR. Otherwise it is the DR parameter at 
the scale Mz- 

4 Note, however, that the conversion is done through a counterterm that involves, as required for the order 

O(atcts) corrections computed in NMSSMCALC order 0(a s ) corrections and no EW corrections. 
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Figure 6: Upper: The partial decay width F(A 2 —>• titi) as a function of at LO (blue/lower dotted), 
including the NLO QCD (green/upper dotted), the NLO EW (red/dashed) and both the EW and QCD corrections 
(black/full). Lower: The relative correction A = (r^ LO — r LO )/T LO in per cent for X =QCD (green/dotted), 
EW (red/dashed) and the full NLO corrections (black/full). The pink line shows the position of the parameter 
point defined in Eqs. (l6371l iKTHTll . 


The leading order width obtained in this scenario amounts to 

T l0 (^2 -> tih) = 15.72 GeV , (6.64) 

where again we have summed over both charge conjugated stop pair final states. The LO width 
differs by 2.7% from the value obtained at tree level with NMSSMCALC, where the Fermi constant 
Gf instead of a is used as input parameter. 

In Fig. [6] (upper) we show the partial decay width T{A 2 —>■ tiQ) at LO, including the EW and 
the QCD corrections, and the NLO width with both the QCD and EW corrections taken into 
account, as a function of m^. which is varied around the parameter point defined in Eqs. (16.571) 
(16.611) with rrif^ « 281 GeVO Note, that stops can still be rather light jTOl. 102 : . T04l[108[ri09j . 
down to about 240 GeV for arbitrary neutralino masses |I02 j and even lower taking into account 
the actual t\ branching ratio )TW| . The figure illustrates the effect of the higher order corrections, 
although the thus obtained parameter configurations are not all in accordance with the applied 
constraints anymore. The lower plot displays the relative corrections 

pNLO _pLO 

A = -, V = QCD, EW, QCD+EW (6.65) 

in per cent. The plots show that both the QCD and the EW corrections are significant and 
come with opposite sign. The QCD correction increase the LO width by ~ 40 — 120% in the 

5 The variation of m,f 1 between 170 and 300 GeV corresponds to a variation of A P s between 1371 and 1721 GeV. 
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investigated range, depending on the value of m^, whereas the EW corrections are almost 
independent of and decrease the cross section by 40%. At = 192 GeV the QCD and 
EW corrections are of same size and cancel each other. The full corrections hence increase the 
width between ~ 0 — 80%, cf. Fig. [6] (lower). And for our parameter point the total correction 
is 


ATQCD+ew^ = 41 % (6.66) 

This plot demonstrates that both the inclusion of the EW and the QCD corrections is required 
in order to properly predict the decay width. 
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Figure 7: Same as Fig. [6j but for the branching ratios. 

In Fig. [7] we show the branching ratios corresponding to the widths of Fig. [HI They have been 
obtained by replacing in NMSSMCALC the corresponding width with our loop corrected width@ 
The branching ratio at LO of our investigated parameter point amounts to 

BR lo (A 2 -> t x t 2 ) = 40.8% . (6.67) 

The net effect of the NLO EW and QCD corrections is an increase of the branching ratio by 

ABR qcd+ew (A 2 -> hh) = 20.8% . (6.68) 

In the plot of Fig. [7] we again vary rn^ around the chosen parameter point, illustrated by the 
pink line in the plot. As can be read off the lower plot the total NLO corrections increase the 

6 In NMSSMCALC the SUSY QCD corrections to the decays into squarks, as derived from [42], are taken into 
account. These include improvements in the decays into sbottoms, which are required in parts of the parameter 
space, that are not relevant for us. Furthermore, we include the EW corrections. We therefore consistently turned 
off the corrections implemented in NMSSMCALC in the decays into squarks and included instead our corrections. 
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branching ratio by up to a bit more than 40% in the shown range. Thus this decay remains 
the most important one also after inclusion of the NLO corrections. 

Figure [ 8 ] finally shows the dependence of the higher order corrections on the gluino mass. 
The EW corrections of course do not depend on nig, while the QCD corrections show a very 
mild dependence on the gluino mass, apart from the region around mg ~ 535 GeV. The kink 
that appears here, arises in the G self-energy at the threshold where = nig + m;. 
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Figure 8: Upper: The partial decay width F(t12 — > trtz) as a function of rrig. The color and line style code is the 
same as in Fig. [6j Lower: The relative correction A = (r5f LO — r LO )/r LO in per cent for A =QCD (green/dotted), 
EW (red/dashed) and the full NLO corrections (black/full). The pink line shows the position of the parameter 
point defined in Eqs. (16.57l) - (16.61l) . 

The size of the higher order corrections sensitively depends on the scenario. Thus we find 
for the scenario defined by 

m UR;C R = m d R ,s R = m Qi , 2 = m hi , 2 = m e R ,iiR = 3 TeV, m- tR = 714.25 GeV, 
ttIq 3 = 1035 GeV, nt^ R = 2776 GeV, m ^ 3 = 2156 GeV, m,f R = 1755 GeV, 

A t = 1246 GeV, A U)C = 1347 GeV, A d ^ b = -1651 GeV, A e ^ r = 769.08 GeV, 

Mi = 460.61 GeV, M 2 = 381.55 GeV, M 3 = 2296 GeV (6.69) 


and 


A = 0.552 , k = 0.030 , A K = -173.51 GeV , fi eS = 446.80 GeV , 

tan /3 = 3.005 , M H ± = 1460 GeV , (6.70) 


resulting in 


Ma 2 = 1461 GeV , = 353.02 GeV and mi 2 = 927.56 GeV (6-71) 
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and an LO decay width of 


(6.72) 
10 .2%, respec- 


T lo (A 2 -> ht 2 ) = 15.24 GeV 

NLO QCD and EW corrections that amount to A^ CD = 23.4% and A EW = — 
tively, resulting in a total correction of 

a QCD+ew^ 2 = 132 % ( 6 .73) 

This can also be inferred from Fig. [9] which shows the NLO corrections to the decay widths as a 
function of m £ . As demonstrated in the lower plot, the NLO QCD corrections are of the order 
of 20-25%, while the EW corrections range between about -17% and -8%, leading to an overall 
increase of the cross section between 3% and 17% due to the combined NLO corrections. 
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Figure 9: Same as Fig. [U] but now for the initial scenario (marked by the pink line in the plot) given by Eqs. (16.691) 
( 10 ) 1 ) 

6.2 Stop Decays into a Pseudoscalar 

Decays of the heavy stop into a pseudoscalar and the light stop, 

t 2 —> Ai + t\ , i = 1,2, (6.74) 

can occur and become important when there is a large mass splitting between the two stop mass 
eigenstates. Stop decays into a SM-like Higgs boson final state have recently been discussed 
in [IIP) and the production of NMSSM Higgs bosons in squark cascade decays in [5511111 j . The 
LO decay width is obtained from Eq. (13.1811 by adapting the kinematic factor and dividing by 
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the color factor 3 and the factor 2 due to the summation of the two charge conjugated stop pair 
final states, 


r LO (t 2 -A- Ait\) — 


A 1/2 (m? ) 


h _ 

I67rmf 

C 2 


11 E Z « G 

3 = 1 


(6.75) 


The NLO SUSY-QCD decay width can be written as 

r NLO (t 2 -A- Aiti) = r LO (f 2 —> Aiti) + TQQ D (f2 —> Ait\) 

where 


^qcd(^ 2 ~^ Ait\) — Re 


with the correction factor 




_ f_2 _ 

247rmf 

C 2 


E Z « G 

J =1 


12 * 1 &s 

A i ■ ^ 

7T 




\k= 1 


a qc D = a^ +a ct A R 

k,t2 k,t2 k,t2 k,t2 


(6.76) 


(6.77) 


(6.78) 


The virtual correction^ and the counterterms are given by the same expressions as for the 
pseudoscalar decay presented in section 01 i.e. 


A v - A v 
A k,i 2 ~ ^ A k 


and A?? = A ? t 


*k,i 2 — A k ■ 

In the real corrections, though, the roles of A^ and f 2 have to be interchanged, 

A R_ = aK 0/2 , . _2 


^ t - 2 = A 


(6.79) 


(6.80) 


with given in Eq. (14.411) . The NLO SUSY-EW corrections are composed of the same 
electroweak virtual corrections to the verted, the same counterterms and the same mixing 
contributions of the pseudoscalar with the Z and Goldstone boson, A^ v and 5 M^r as in the 
pseudoscalar decay. In the real corrections, however, again M\. and m ~ have to be interchanged. 
We hence have 


A 1 / 2 


AA) ~ 1 form? 


t2 


E Z « G 


3 = 1 


+ E z « G i 2 ,* E Z <*Y 


h 


\3 = l 


\k=1 


-2Re E TAiPX E + SM > 


G,Z 

mix,i 


0=1 


\k=1 


with 


A 1/2 = A 1/2 (m|,M|.,m?) 




(6.81) 


(6.82) 


7 The diagrams are the same as in Fig. [I] but with fe in the initial and Ai in the final state. 

8 In the Feynman diagrams of Fig. [3] simply the t 2 leg has to be crossed to the initial and the Ai leg to the final 


state. 
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Figure 10: The decay width V(t 2 —> A\ti) as a function of at LO (blue/upper dotted), including the NLO QCD 
(green/lower dotted), the NLO EW (red/dashed) and both the EW and QCD corrections (black/full). Lower: 
The relative correction A = (Tx LO — r LO )/r LO in per cent for X =QCD (green/dotted), EW (red/dashed) 
and the full NLO corrections (black/full). The pink line shows the position of the parameter point defined in 
Eqs. (1051) - (1061) . 


and 


R, EW 
^ k,t 2 


167r 2 


A % 


(Ml Om|) 


(6.83) 


and A^ v and S M^r l given in section [5j The full NLO decay width including SUSY-QCD and 
-EW corrections is given by 


pNLO, full _ pNLO I p(l) / f? o r\ 

1 - 1 EW + 1 QCD i 

with rQp D given by Eqs. (16.771) (16.801) and by Eq. (I6.81|) . 

In order to show the impact of the SUSY-EW and -QCD corrections on the stop decay width 
we chose the following parameter set, which leads to an NMSSM Higgs and SUSY spectrum in 
accordance with the LHC data: 


m UR,c R = m d R ,~s R = m Qi,2 = m Li,2 = = 3 TeV > m t R = 748 - 07 GeV » 

rriQ^ = 1259 GeV, = 1709 GeV, = 1637 GeV, rrif R = 1618 GeV, 

A t = 1589 GeV, A u , c = 1675 GeV, A dtStb = -669.04 GeV, A e ^ T = 179.93 GeV, 

Mi = 645.51 GeV, M 2 = 272.11 GeV, M 3 = 2511 GeV , (6.85) 

and the NMSSM specific input parameters 

A = 0.588 , k = 0.378 , A K = -675.73 GeV , /r e ff = 385.90 GeV , 

tan J3 = 1.529 , M H ± = 639.08 GeV . (6.86) 
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This results in the two-loop corrected mass Ma 1 of the lighter pseudoscalar A± and the stop 
masses m^ , 

M / \ l = 637.40 GeV , rn^ = 342.76 GeV and mi 2 = 1153 GeV . (6.87) 

And for the stop soft SUSY breaking trilinear coupling in the OS scheme we get 

A° s = 1675 GeV . (6.88) 

Figure [TO] (upper) shows the decay width r(Q — > A\t\) at LO, including the NLO EW, the 
NLO QCD and both NLO corrections, as a function of mi , varied around the chosen parameter 
point, marked by the pink line in the plots. The LO decay width reaches 

r LO (f ~ 2 -> AiU) = 9.26 GeV . (6.89) 

In the whole investigated range, the EW and QCD corrections are significant, decreasing 
together the LO width by ~ 12 — 15%. Both corrections are of similar size, where, depending 
on the parameter point, once the QCD, once the EW corrections are more important. At our 
starting parameter point both corrections are almost equal resulting in a decrease of 

ATQCD+ew = _ n 7% ^ (6.90) 

cf. Fig. [TO] (lower). Our results demonstrate that also in the stop decays both the QCD and EW 
corrections have to be considered for a meaningful prediction of the stop decay width. 

7 Conclusions 

The search for New Physics is one of the main tasks at the LHC. In the absence of any direct sign 
of new resonances so far, the precise investigation of the Higgs sector becomes more and more 
important. Physics beyond the SM might reveal itself in modified Higgs decay rates compared 
to the SM expectations or in the discovery of additional Higgs bosons, unambiguous sign of a 
non-SM Higgs sector. In view of the complexity of the experimental analyses and the plethora 
of still possible New Physics extensions, it is evident that the success of this research program 
depends on the precise predictions of parameters and observables from the theory side. In 
this paper we have calculated the NLO SUSY-EW and -QCD corrections to the decays of a 
pseudoscalar NMSSM Higgs boson into stop pairs and of the heavier stop into the lighter stop 
and a pseudoscalar Higgs boson. Both processes can become important in certain regions of the 
parameter space and hence contribute to the discovery channels of either the pseudoscalar Higgs 
boson and/or the stop quarks. The NLO corrections turn out to be important and, depending 
on the scenario, the EW corrections can be of same size as the QCD corrections and also come 
with opposite sign. Therefore not only the inclusion of higher order corrections is important 
but also the consideration of both the QCD and the EW corrections is indispensable for making 
reliable predictions. 
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Appendix 

A The Loop Functions 


The D dimensional one-loop integrals encountered in the calculation are the scalar one-, two- 
and three-point functions Aq, Bq and Co as well as the coefficient of the two-point tensor integral 
of rank one, B\. They are defined as 


A 0 (m) 


B 0 (p 2 ;mi,m 2 ) 


Co{pl,pl,pl 2 ;rni,m 2 ,m 3 ) 


p^Bi{p 2 ]mi,m 2 ) 


where 


16ttV _D J 
16vr 2 p 4 ~ D J 
167tV 4 - D J 


dP q 1 

i(2n) D ( q 2 — m 2 ) 
d D q _ 1 

i(2ir) D (q 2 — m\)[(q + p) 2 
dPq 

i(2ir) D 


1 



( q 2 - mf)[(q + pi) 2 - mj \[(<? + p l2 ) 2 - m§] 


16vrV 4 “ D 


d D q 


% 


i(2ir) D ( q 2 — m\)[(q + p) 2 — m 2 ] 


(A.l) 

(A.2) 


(A.3) 

(A.4) 


P12 = Pi + P2 


(A.5) 


with p,pi,p 2 denoting the external momenta, which are taken as incoming, and m,mi 1 m 2 ,m 3 
the masses of the loop particles. 
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